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The properties of the quark and hadron Universe are explored. Kinetic 
theory considerations are presented proving that hadron abundances after 
phase transformation from quarks to hadrons remain intact till abundances 
of hadrons become irrelevant. The hadronization process and the evolution 
of hadron yields is described in detail. 

1. Motivation: The big-bang and laboratory experiments 

1.1. A short history of quark and hadron Universe 

At early times (5 < t < 20 //s) in the standard model of the big 
bang, matter as we know it was dissolved in a state known as the quark- 
gluon plasma (QGP), in which leptons, gauge mesons, and quarks u,d,s 
were deconfined and propagated freely — all hadrons are dissolved into their 
constituents. At that time, strongly interacting nearly massless particles 
u, d, s, G controlled the fate of the Universe: allowing for their interaction 
these degrees of freedom comprised nearly 70% of all energy and pressure. 
The remaining 30% was shared by f ejAljT , e, /i, 7 pp. 

The epoch we address begins just before the quark-gluon Universe turns 
into the hadron Universe and ends when we reach the era of the electron, 
positron, neutrino, and photon dominated Universe. Our interest is to de- 
scribe in detail how hadronization in the Universe happened and to explore 
the Universe's properties while it expands towards neutrino decoupling. In 



* Presented by JR at 52 Krakow School of Theoretical Physics: Astroparticle Physics 
in the LHC Era, Zakopane, May 19-27, 2012 



(1) 



2 



Rafelski'QGU3 printed on November 20, 2012 



outline, this epoch proceeds as follows: Around the hadronization (Hage- 
dorn) temperature T/j ~ 160 MeV, a phase transformation to the confining 
state occurs, and matter as we know it today forms. In the first moments 
this hot hadronic 'gas' phase (HG) is nearly symmetric in its matter and 
antimatter components, and it still comprises about 70% of all energy, so 
that for a short instant the Universe is hadron matter dominated. It takes 
about 0.1 s for hadrons and antihadrons to annihilate so thoroughly that 
in essence only the net baryon number we see around us today remains: 
baryons in the present day Universe retain a residual 4.6% energy contri- 
bution, but only 10~ 10 fraction of particle number. Therefore, when the 
Universe was at a temperature prior to nucleo-syntheis but after neutrino 
freeze-out, let us say near T = 2 MeV, this residual baryon matter fraction 
was not larger than 10 -7 and the Universe was dominated by photon and 
neutrino 'radiation'. 

There is consensus that at approximately 1 ms (T ~ 35 MeV) , the 
matter-antimatter annihilation was nearly complete. What we do not see in 
the literature is that the truly interesting feature of the early Universe in the 
following epoch is an enormous remaining pion component, with an abun- 
dance 10 7 exceeding that of nucleons, and which becomes subdominant to 
nucleons only at T = 5 MeV. Similarly, strangeness is present and becomes 
subdominant in abundance to baryons only at T = 20 MeV. Strangeness, 
and in particular CP asymmetric neutral kaons, are recreated at a time 
scale which is hadronic while the Universe evolves. There are many reac- 
tions potentially enhancing the outcome of a small CP breaking in the kaon 
system, prompting the question, could we have too quickly discounted the 
possibility of baryon asymmetry originating in the hadron Universe? To 
answer this question we must first find an accurate and precise description 
of the Universe in this time period. 

This report presents a brief synthesis of our effort to improve the un- 
derstanding of the early quark and hadron Universe which we pursued over 
the past decade. Our efforts were initiated in the study of time scales of 
the transformation from QGP to HG phase [2]. We soon added to this 
a more detailed description of the Universe properties [3], assuming that 
chemical equilibrium prevails. Later, using kinetic theory we showed that 
the expansion of the Universe is indeed slow enough to assure that chemical 
equilibrium prevails [H [5] . 

1.2. Laboratory experiments addressing the early Universe 

Relativistic heavy ion collision experiments at RHIC and the LHC pro- 
vide an important opportunity to simulate some aspects of these early Uni- 
verse conditions. In particular, in most central heavy ion collisions, a small 
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nuclear-sized volume attains energy density and temperature comparable to 
those that prevailed in the QGP in the early Universe, this is the micro-bang. 
The particles observed later in the detectors are products of the following 
rapid, explosive hadronization, an energy-to-matter conversion analogous 
to the one occurring very slowly in the early Universe. Moreover, the QGP 
formed this way in the laboratory has several differences: the characteristic 
timescale of the heavy ion collision is much shorter r ~ 10 -23 s which means 
that only strongly interacting particles equilibrate, and the baryon density 
is much higher, though baryon density attained at LHC is now so low that 
we have yet to measure it. Taking into account these differences, we show, 
in this report, how the tools developed to study 'micro-bangs' are applied 
to the early Universe hadronization and subsequent annihilation dynamics. 

An important question we can address quantitatively is, what conditions 
in the early Universe QGP phase yield the observed matter asymmetry? The 
present-day small value of the baryon-to-photon ratio (which quantifies the 
amount of residual matter - see section 13. 1[) is seen by many as the result 
of near-complete annihilation following hadronization. This is supported 
by the homogeneity of the Universe which is seen as being consistent with 
absence of antimatter on large scales, and some workers think that separa- 
tion into matter-antimatter domains on a scale smaller than the observable 
Universe is unlikely [6j [7] . 

However, more recent views on the subject sound different [8j, we cite: 
"Though observation strongly restricts possible existence of antimatter do- 
mains and objects, . . . , it is still not excluded that antimatter may be 
abundant in the Universe and even in the Galaxy, not far from us. That is 
why there is an active search for cosmic antimatter . . . " . A prominent ex- 
ample of active antimatter search is the AMS experiment [9] running at the 
International Space Station, searching for anti-a. Discovery of anti-a would 
provide the first evidence that matter-antimatter separation has occurred. 
Such separation could arise in the early Universe in slow hadronization of 
QGP. We study this process and quantify the electrical charge distillation 
in Section 13.31 

We proceed starting in Section [2] by first recalling the statistical physics 
of relativistic quantum particles, and how it is used to describe the prop- 
erties and transition from the QGP to HG in the early Universe. Special 
attention is given to the understanding of chemical potentials which play 
a pivotal role in the understanding of particle abundances: in a Universe 
in thermal and chemical equilibrium the matter-antimatter asymmetry is 
expressed by non-zero values of the chemical potentials. We study their 
magnitude quantitatively in Section [3] assuming that the Universe is on av- 
erage charge neutral, the lepton number is equal to the baryon number, and 
entropy per baryon is derived from baryon-to-photon ratio. 
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Our study of particle abundances relies on hadron reactions occurring 
even in relatively low T environment at sufficient speed to assure the pres- 
ence of all hadrons and we establish this considering key relaxation times 
in section SJ The key reaction are introduced, which connect hadrons with 
highly abundant photons and leptons, and equilibrate with the surviving 
pions. The final Section [5] concludes. 

2. Statistical hadronization and the early Universe 



A key feature of hadronic and QCD interactions is that the interac- 
tions are strong, so that at moderate temperatures or densities collisions 
are frequent. In fact, they are frequent enough that statistical equilibrium 
is (partially) attained even within the duration of relativistic heavy ion col- 
lisions At ~ R/c = 5-8 fm/c. The temperature and energy density achieved 
in the collision are sufficiently high that quarks and gluons are expected to 
be 'deconfined' and the relevant degrees of freedom across the volume in- 
volved in the collision; in other words, a thermalized drop of QGP is created 
in the lab. This dense and highly compressed state of matter expands in 
a 'micro-bang' and as temperature drops, hadronizes. The produced yield 
of hadrons relates is described within a statistical hadronization approach 
which relies on use of phase space size to predict particle abundances. 

In consequence, relativistic heavy ion collisions evolve via several steps, 
summarized here for comparison to the evolution of the early Universe 
plasma [2]. 

1. Initially, a small thermalized quark-gluon 'fireball' is formed from the 
collision. Thermalization is very rapid in the QGP; the timescale for a 
species i is approximately its mean scattering time with all other species j 



where <Xjj is the cross-section for energy-exchanging interactions, Vij — > c is 
the relative velocity of massless components and rij is number density. For 
typical QGP parameters, 

Ttherm.i = 0.2 — 2 fm/c, rij ~ 2 — 10 fin ^ , aij ~ 2 — 5 mb. (2) 

This estimate agrees with the thermalization timescale r t h e rm ^ 10~ 23 s 
determined from particle spectra and yields. The microscopic mechanisms 
leading to such rapid thermalization are still the subject of intense study. 
Even so, we expect this result to be valid qualitatively in the early Universe, 



2.1. Statistical hadronization, qualitatively 
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and, comparing to the early Universe timescales below, implies the primor- 
dial QGP is thermalized. 

2. Chemical equilibration of quarks u, d and later s takes somewhat longer, 
~ 1.5 fm/c for u, d and ~ 5 fm/c for s. In general, chemical equilibra- 
tion timescales are longer than thermalization timescales, due to smaller 
cross sections for particle (quark pair) production as compared to energy- 
exchange. However, detailed studies of the QGP properties suggest that 
time is sufficient for chemical equilibration of u, d, s at LHC. However, the 
heavy quarks c, b do not have time to come to chemical equilibrium due to 
their higher masses implying longer chemical relaxation times. Remarkably, 
their abundance at in LHC heavy ion collisions is expected to be well above 
chemical equilibrium due to production in initial high energy parton col- 
lisions. The same kinetic theory approach that produced relaxation times 
for strangeness, also allows us to compute the relaxation time of c, b quarks 
and we find that heavy quarks are in chemical equilibrium at the time scale 
(microseconds) governing the end phase of QGP in the Universe near to 
T ~ 160 MeV. 

3. Hadronization occurs as the fireball expands and its temperature drops. 
In the QGP hadronization that follows on its formation at LHC, the matter 
expands into empty space explosively and thus the hadronization process is 
too rapid to allow re-equilibration of the final state hadrons. When chem- 
ical equilibrium of a particular flavor can no longer be maintained in the 
expanding fireball, the flavor 'freezes out' — its particle number in the co- 
moving volume remains approximately constant while its density is diluted 
by the expanding the volume. Here, the early Universe differs; hadroniza- 
tion encompasses all the volume and is expected to be much slower with 
chemical equilibrium maintained throughout. 

The dynamics of chemical equilibration is especially important in rela- 
tivistic plasma, in which energy can be converted to and from particles and 
antiparticles. Note that there are two distinct senses of chemical equilib- 
rium: 

i) Absolute chemical equilibrium is the level to which energy is shared into 
accessible degrees of freedom. A black body photon gas achieves absolute 
chemical equilibrium almost instantaneously because the photon is massless 
and photons are equilibrated easily at the hot walls of the oven. However, a 
massive particle whose particle number is initially zero takes time to achieve 
chemical equilibrium due to the necessity of finding particles able to react 
and create the massive particle. 

ii) Relative chemical equilibrium reflects the distribution of an already exis- 
tent element or component among different compounds. Relative chemical 
equilibrium is thus associated with the 'usual' sense of chemical potential 
ensuring the conservation of conserved 'charges', such as baryon number. 
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In order to be able to describe in time the approach to absolute chemical 
equilibrium it is necessary to introduce, aside of chemical potentials such as 
Hb for baryons, additional abundance fugacity 7/(i) controlling the absolute 
yields of pairs of particle / so that the total fugacity, e.g., for nucleons and 
antinucleons is 

T N = lN (t)e^ T , T^ = lN (t)e^ T , Wb = -^b. (3) 

Note that a quark is always produced paired with an antiquark so that the 
net baryon number is not changed. The change in chemical potential due to 
an antiparticle is for this reason opposite in sign to that due to a particle as 
indicated above; provided, of course, the overall yield is scaled appropriately 
by the fugacity j(t) taking into account the approach to absolute chemical 
equilibrium. 

The consequences for the energy budget are seen in the first law of 
thermodynamics 

dE = - PdV + TdS + TlnT f dN + TlnT w dN, (4) 
= - P dV + T dS + n b (dN - dN) + Tlnj N (dN + dN). (5) 

While in pair production the coefficient of \x\, vanishes, the energy required 
to change the number of nucleon-antinucleon pairs is related to 7/v- Near 
to absolute chemical equilibrium 7^ — >■ 1, the last term vanishes; small 
fluctuations in the number of pairs do not influence the energy balance. 

Statistical hadronization incorporates these dynamics of the quark and 
hadron 'chemistry' to provide quantitative predictions of the statistical 
properties of the final state HG. It successfully describes particle produc- 
tion in heavy ion collisions [2], and, as a general framework derived from 
the underlying theory of QCD and phenomenology of the HG phase, we can 
apply it to study the early Universe hadronization and the following HG 
dynamics. 

2.2. Statistical models of QGP and HG 

Thermodynamic properties of the QGP and HG phases are computed 
from the partition functions ln^QCP and In Zhq; f° r details see [2], we 
give here a short introduction. An important difference to usual class room 
circumstance is that particle numbers vary. This possibility of evolving 
changes in (quasi-)particle number is accounted for in the grand canonical 
partition function Z 

Z(p,V,T f ) =vJ29fj-0f (ln(l±T/6-*/) +ln(l± !>-*/)) , 

(6) 
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which is a function of the volume V, inverse temperature (3 and the fugacity 
Tj j which accounts in its properties Eq. ([3]) for the possibility that absolute 
chemical equilibrium is not reached. The total partition function for each 
component is sum over the (quasi-)particles / present. Here, e 2 = p 2 + m 2 
is the energy. 

For the QGP in absolute chemical equilibrium the momentum integral 
can be carried out when the mass is 'small' 

|ln^ QGP =-S + i |^ Cl (vrr) 4 



t w 157r2 

j=u,a,s 



15a s 50a s 2a s 
ci = 1 h • • • , c 2 = 1 - — h • • • , c 3 = l h • • • , 

47T 217T 7T 

which describes quantum gases of quarks and gluons, including also the first 
order perturbative QCD corrections, see the a s terms in the Cj coefficients, 
and a confining vacuum energy-pressure component B ~ 0.19 GeV fm -3 . 
The temperature provides the scale for the evaluation of ct s (fj,) and is esti- 
mated as fi = 2itT or with finite chemical potential fi = 2^Jir 2 T 2 + fjij. 

The HG partition function, used in our computation, included a sum 
of partial gas contributions from all hadrons having mass less than 2 GeV, 
and we apply finite volume corrections [10]. Together with the QGP-liquid 
model [UJ, this framework is a phenomenological description of QGP equa- 
tions of state which agrees with properties of quantum chromodynamics 
(QCD) at finite temperature obtained in lattice QCD in the limit of van- 
ishing particle density (^tj — > 0) [121 E2J E] • 

Recall that the statistical partition functions imply the thermodynami- 
cal relations. Therefore, the entropy and number density in each phase can 
be obtained by introducing the free energy 

^ = -^mz = -p, (8) 

where P is the pressure. We will require below the entropy density s and 
the baryon number density ng, being 1/3 the quark density 

ldJ 1 1 d T 



See Chapters 4, 10 & 16 of [2\ for more details. 
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2.3. Early Universe conditions 

First, let us check that local thermal equilibrium prevails throughout 
the period of interest using Eq. (fl]). In the HG phase, where n ~ T~ 3 
and particle velocity is non-relativistic v ~ T 1 / 2 , r t h e rm ~ T~ 7 / 2 . Scaling 
^therm ~ 10 -23 s from hadronization means Ttherm ^ 10~ 14 s at T = 1 MeV. 
Comparing to the age of the Universe t ~ 1 s shows that the assumption of 
local thermal equilibrium has a large margin of error for the whole process 
of hadronization and equilibration discussed below. 

In [2], it is estimated that the phase transformation from QGP to HG 
takes Th ~ 10 [is in the early Universe, many orders of magnitude longer 
than chemical equilibrium timescales. Therefore, chemical equilibrium of 
hadrons made of u, d, s quarks is firmly established at the end of the phase 
transformation. Baryon annihilation reactions cease near T = 40 MeV. 
Hadron abundance evolution in the early Universe and possible deviations 
from the local equilibrium will be studied down to T ~ 2MeV in Sec. HI Our 
detailed study here demonstrates that chemical equilibrium prevails down 
to neutrino decoupling. 

An essential characteristic of the expansion and cooling of the Universe is 
that the evolution is isentropic. Combining the first law of thermodynamics 
Eq. ([U with the Friedmann equation, there is a direct relation between the 
total energy density of the system and the scale factor of the Universe 

(t is the comoving coordinate time in Friedmann-Lemaitre-Robertson- Walker 
coordinates). For this reason, the scale factor of the Universe will remain 
implicit, and the system (=the Universe) will be tracked in terms of its 
thermodynamical properties T,fif. T is effectively the clock tracking the 
progression of the system from QGP through hadronization to HG and 
subsequent equilibration, and, in the next section, a system of constraints 
relating the [if will be developed and solved at each T. The explicit value of 
the expansion rate dR/dt is needed only later in Sec. 0J when reaction rates 
must be compared to dR/dt to determine freeze-out times/temperatures. 

The framework so far is fixed by well-established conditions of the early 
Universe expansion and the statistical models of the QGP and HG phases, 
which are motivated by experiment. However, we should mention several 
assumptions implicit in this setup and relevant to cosmology: i) dark energy 
is irrelevant at this stage of the Universe's evolution; ii) dark matter does 
not affect the populations of visible (Standard Model) particles, neither by 
its decay /freeze-out nor by secondary effects; and iii) there are only left- 
handed neutrinos and corresponding antineutrinos. 
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3. From QGP to cold hadron Gas: 200 > T > 2 MeV 

In a system of non-interacting particles, the chemical potential [if of 
each species / is independent of the chemical potentials of other species, 
resulting in a large number of free parameters. In the early Universe QGP 
and HG, the number of independent parameters is reduced by the many 
chemical particle interactions occurring between the species. 

First, in thermal equilibrium, photons assume the Planck distribution, 
implying a zero photon chemical potential: /x 7 = 0. Second, for any reaction 
VfAf = 0, where Vf are the reaction equation coefficients of the chemical 
species Af, chemical equilibrium occurs when Vf[if = 0, which follows from 
minimizing the Gibbs free energy. 

Next, chemical and thermal equilibrium means that particle-antiparticle 
pairs can be created from and annihilate into the thermal bath of photons, 
i.e., the reaction / + / ^ 2j proceeds freely in both directions. There- 
fore, [if = —fJ>f whenever chemical and thermal equilibrium is attained. In 
particular, when the system is chemically equilibrated with respect to weak 
interactions, the following relationships hold, see, e.g., [IB] : 

Me - Hu e = [i-e - 1% = Mr ~ fh/ T = A/4, (11) 

A i « = Md-A^, n s = n d . (12) 
The baryon chemical potential is 

3 3 

In relative chemical equilibrium, the chemical potential of hadrons is equal 
to the sum of the chemical potentials of their constituent quarks. For ex- 
ample, Y,°(uds) has chemical potential /i s o = [i u + [id + [i s = 3 [id — A/i/. 

Finally, neutrino oscillations, especially with the large mixing angle now 
strongly supported by experiment [16j . imply that neutrino number is freely 
exchanged between flavors v e ^ ^ v T and hence 

[iu e = ih,,, = /V = Hv (14) 

Together, these conditions reduce the number of independent chemical 
potentials to three. We choose to track fid, He, and [i v . 

3.1. Three constraints 

When the Universe is in a single phase (QGP or HG phase), the three 
free chemical potentials are determined by the following three criteria: 
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i. Charge neutrality {Q = 0) is required to eliminate Coulomb energy. 
This is written as 

/ 

where Qi is the charge of species /, and the sum is over all particle 
species present in the considered particle phase. 

ii. Net lepton number equals net baryon number (L = B) is phenomeno- 
logically motivated in the context of baryogenesis. This leads to the 
constraint 

n L -n B = ^2(L f -B f )n f (fi f ,T) = 0, (16) 
/ 

where Lf and Bf are the lepton and baryon numbers of species /. This 
condition is of course not proved by experiment and future studies 
must consider drastically different variants. 

iii. Constant entropy-per-baryon (S/B) is equivalent to the statement that 
the Universe evolves adiabatically, and hence can be written 

s _ Ef Sf(t*f,T) 
— = -r = const, (17) 

where sj is the entropy density of species /. 
The constant s/ns is estimated from the value of the ratio 

v = —, mo = io lc V (is) 

n B 

At 100 > T > 2 MeV, there are 43/4 degrees of freedom [1] due to 
photons, e, /i, r-neutrinos, and electron-positron pairs, which provide the 
dominant fraction of the entropy and pressure, controlling the dynamics of 
the Universe and thus the speed of expansion at the time big bang nucle- 
osynthesis (BBN) occurs. We are interested in the Universe entropy content 
per baryon 

s 1 v / s 7 n v s v n e s e \ n 7 
— = — V **=(— + + — , 19 

nB riB n^n v n y n e J ub 

In the Standard Model, with only left-handed neutrinos, the relative num- 
ber of degrees of freedom of each species is known and relates the number 
densities by (4/3)n e = 2n 7 and (4/3)n;, = 3n 7 . The factor 4/3 is the bose 
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Fig. 1. The fit of rj and A7V„ from BBN (deuterium and 4 He abundances) obtained 
in Figure 4 of [T7]. The major axes of the 68% (solid) and 95% (dashed) confidence 
ellipses align with the overlaid contours of constant s/tib- Each contour means a 
change of s/hb by 0.5 10 10 . 

to fermi ratio of particle densities arising from the quantum statistics of the 
massless quanta. The entropy per particle for a boson is (s/n)boson = 3.601 
and for a fermion is (s/ n )fermion = 4.202. Inserting this into Eq. (fl~9|) we find 
the known outcome. 

However, recent BBN analysis combined with the cosmic microwave 
background (CMB) fluctuations have suggested that the number of rela- 
tivistic degrees of freedom present at T > 1 MeV is somewhat higher than 
the Standard Model expectation (T7J. This possibility is discussed in terms 
of extra neutrino degrees of freedom 3„ — > 3 U + AN U , the AN u -r/ fits from 
BBN, see circles (red) showing the best fits in Fig. [TJ taken from Figure 
4 of LIT] - The BBN best fit values [17] arise for AN U = 0.66±g;|| and for 
v = n B /n^ = 6.27 ± 0.34 10~ 10 . 

We generalize Eq. (|19|) to variable neutrino abundance and obtain 




(20) 



(19.35 + 3.15AAU^. 

n B 



(21) 
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Note that this relation applies at/after BBN. We assume that there is no 
change in entropy content thus s/ub we find at time of BBN applies in the 
Universe prior to BBN. Baryon number is conserved. 

We are interested to understand what this means for the variable s/tib 
which enters our evaluation of the Universe properties. In Fig. [IJ we 
overlay contours of constant s/tib on Figure 4 of |17| with lines of con- 
stant s/riB- We recognize that the relevant thermodynamic variable in 
these considerations is the entropy per baryon, and this analysis leads to 
s/tib = 3.42^Q3g 10 10 . In the following figures, however, the numerical re- 
sults shown are based on an older value obtained for s/ng = 4.5^}'^ 10 10 
based on similar analysis of 10 years ago [18] . Since this value of s/ub is 
consistent within its la expectation with the updated 2012 value, we can 
expect the quantitative outcome changes only slightly. 

3.2. Chemical potentials and particle abundances in early Universe 

For each temperature T, the conditions Eqs. ()15p - (jl7p form a system of 
three coupled, nonlinear equations of the three chosen unknowns (fid, l^e, 
and fi u ). These equations were solved numerically [3] using the Levenberg- 
Marquardt method [19] and the results shown in Fig. [2j The bottom axis 
shows the age of the Universe and the top axis the corresponding tempera- 
ture. At low temperature, ^ approaches (weighted) one-third the nucleon 
mass (2m n — m p )/3 = 313.6 MeV reflecting the dominance of protons and 
neutrons in their classical Boltzmann limit T <C mj. 

From Fig. [21 the value of the baryon chemical potential just before the 
phase transition is 

fi B = 0.33l°;^ eV. (22) 

More generally, we can say with near certainty that the values of the chem- 
ical potentials required to generate the observed matter-antimatter asym- 
metry are 

Hf <0.5eV, f = d,e,v. (23) 

This statement admits a generous margin of error for uncertainty in the 
transition temperature 2), ~ 160 MeV as well as the displayed error bars, 
which show the effect of the uncertainty in s/ub obtained from r]. These 
error bars could be improved using the updated measurements of rj [20] . 

From the solution for the chemical potentials, the evolution of particle 
numbers follows. These are shown in Fig. O Strangeness, present in kaons, 
persists in significant number down to T = 10 MeV. Pions are even more 
abundant, their density falling below density of nucleons only at T ~ 6 MeV. 
Among the leptons, note that muons also persist to relatively low tempera- 
tures in the HG phase. Pion and muon evolution and chemical equilibration 
are discussed in more detail in the next Section [JJ 
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Fig. 2. Chemical potentials fid, He, and fi v tracked from the QGP through the 
phase transformation to the HG. The error bars correspond to the uncertainty in 
n, see below Eq. (fT?)) . At right, the time immediately around the transformation is 
zoomed. 

Fig. H] summarizes by showing the fraction of the luminous energy (den- 
sity) in baryons as a function of temperature. At the QGP-HG phase trans- 
formation, the fraction of energy in baryons is ~ 10%, a local maximum 
before annihilation rapidly reduces it. The baryons contribute significantly 
to the energy budget again only once the Universe has cooled and entered 
the matter (rest mass) dominated era. 

3.3. Mixed phase in QGP-HG transformation 

It is by now well established that at low baryon density the QGP to HG 
phase change is a smooth transformation. The macroscopic conditions i-iii. 
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Fig. 3. (Color online) Top: The evolution of hadron densities after hadronization. 
(total) indicates the sum of K+ , K- , K° and K . Note pion-nucleon 'equality' at 
T ~ 6 MeV. Bottom: The evolution of lepton densities after hadronization, and 
(tot.) indicates the sum over neutrino flavors. 
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Fig. 4. The fraction of the luminous energy in light particles (photons and leptons) 
mesons and baryons. Note that temperature decreases from right to left. 



above are satisfied on average and for the system as a whole, but can be 
violated locally. 

To solve the conditions Eqs. (|15p - (|17p in the simultaneous presence of 
both phases, the total partition function is parametrized as 



lnZ tot = /hg In^HG + (1 - /hg) IiiZqgp, 



(24) 



in which /hg represents the fraction of total phase space occupied by the 
HG phase. For example, Eq. (|15p is now generalized as 



Q = = ng GP Vqgp + Vhg 



HG 



(25) 



tot 



(1 - / HG ) ™§ gp + /hg 



n, 



HG 

Q 



and analogous expressions hold for Eqs. Q16|) and (|17p . The total volume 
Vtot is n °t relevant to the solution. 

Solving Eq. (|25p and its companions determines baryon number fraction 
in each phase as a function of /hg> shown in Fig. [5j The QGP preserves a 
larger fraction of the baryon number density throughout the transformation, 
and the ratio is nearly constant n^ GP /n^ G ~ 3. 

With the additional assumptions that /hg evolves linearly in time and 
that the total duration of the phase transformation is 10 fis, we obtain the 
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f HG 

Fig. 5. The fraction of baryons in each phase as a function of the parameter /hg, 
which is the fraction of total phase space occupied by the HG phase. 

solution for the Universe during the transformation, indicated by the vertical 
lines in Fig. In reality, both assumptions are sensitive to the equations of 
state of each phase and the dynamics of the phase transformation. A more 
complete solution would be obtained from transport theory applied to the 
coexisting phases. 

From this solution, the net charge per baryon uq/hb is calculated in 
each phase as a function /hg> which is independent of the additional as- 
sumptions. Protons and neutrons being the lowest excitations in the HG 
phase, the HG takes on a positive charge as soon as the transformation 
begins. The QGP therefore takes on a negative charge density, which is ini- 
tially tiny since it occupies the larger volume, yet it can cause large variation 
in local electric potentials. This charge distillation arising as a dynamical 
asymmetry has been discussed in the context of strangeness separation and 
strangelet formation |21|, [22] 

The sign of the charge distillation remains the same throughout the 
transformation, and therefore the total charge of the remaining QGP is 
increasingly negative as the transformation proceeds. One expects the re- 
sulting electromagnetic potential to alter the chemical potentials for charged 
species and thereby to have a feedback effect. Flows of charged particles 
will alter the uniformly small net baryon density and thus any local initial 
baryon-antibaryon asymmetry. Resolving these dynamics clearly requires 
transport theory, as an example of the remark above. 
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Fig. 6. Net charge (including leptons) per net baryon number in the HG and QGP 
as function of the HG fraction /hg- Horizontal line at zero guides the eye. 



4. Pion and lepton equilibration: 50 > T > 2 MeV 

We now move to a closer study of the period of annihilation, decay and 
chemical equilibration immediately following the QGP-HG phase transfor- 
mation. As shown in Fig. 01 this involves a large number of mesons produced 
during hadronization which subsequently decay, affecting the lepton popu- 
lations. The point we address here is that the same reaction that drives 
decay, by detailed balance, also recreates the decaying particle. Here of 
great significance is that the Universe is filled with a background gas of 
photons and leptons that act as a buffer to all hadron decay processes. 

We demonstrate the importance of the background gas considering the 
pivotal reaction 

vr° f± 7 + 7, (26) 

which keeps hadrons in chemical equilibrium throughout the Universe evolu- 
tion. This outcome is in contrast with the naive expectation that tt° disap- 
pears, based on considering only the 7r° lifetime in vacuum tq = 8.4 10~ 17 s 
to the much longer time scale of the Hubble rate of expansion 

H - 1M ^r (omiByl (^mev) 2 - (27) 

in which g* is the number of degrees of freedom contributing to entropy 
at the time (fermionic contributions are modified by the factor 7/8) and 
Afpi = 1.22 10 19 GeV is the Planck mass. 
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Note that the large pion abundance that goes along with reaction Eq. (|26p 
can influence the relic neutrino density and potentially even BBN which im- 
mediately follows this 'pion' era - we hope to return to these topics in the 
near future. Reaction Eq. (|26p is example of important family of one-to-two 
reactions which have only be studied recently [U [5] . Another example is the 
decay of charged pions ir^ can impact the residual populations of muons, 
electrons and neutrinos both through decay and scattering 

7r± - /i± + v^vj. (28) 

The decay of muons also contributes to neutrino density 

M ± ^ e ± + v e (y e ) + v^). (29) 

Two-to-two reactions maintain equilibrium within and between charged 
pion populations 

7T + vr — 7T+ + 7T~, 7 + 7 — 7T + + 7T~, (30) 

and lepton populations 

e + + e -^ M + + M - ) 7 + 7 - i+ + r, 

ir + + n- ^/+ + r, (l = H,e). 

In the following we describe the chemical relaxation times of each of these 
reactions and compare it to the other relevant timescales, especially the 
expansion of the Universe. 

Consider the first example, Eq. (j26|) controlling ir° chemical equilibra- 
tion with the background of thermal photons. The population equation 
describing the evolution of ir° number can be cast into the form a differen- 
tial equation for the 7r° fugacity T3: 

^T 3 = — T 3 + — T 3 + -(T 1 T 2 -T 3 ), (32) 
at t t t s r 3 



see discussion in [24J, partially recounted in the appendix. We keep numer- 
ical subscripts here rather than simplifying with 1 = 2 = 7 so that later 
generalization is clear. Here, the tt and t$ are kinematic relaxation times 
for the temperature and entropy: 

1 3 / (in 3 \ 1 d ( n 3 \ dT 

V T = ~ 9 * T [dTJ dr [t^J W (33) 

1_ _ n, (dr^Y 1 dln(g*T 3 V) dT 
T S ' T 3 UtJ dT df 1 1 



Rafelski'QGUS printed on November 20, 2012 



19 



The leading minus signs mean that tt,t$ > 0. 

T3 is the chemical relaxation time incorporating effects of the background 
gas [5]: 

T3 -vd^ T3 {-lvdT) -vdr-J lT2 [-lvdT) ' (35) 

where N3 is pion number. dW/dVdt is the rate per unit volume for the 
process in subscript accounting for the presence of the background gas of 
reactants and products, and which therefore contains dependence on the cor- 
responding fugacities of the reaction participants. Thus, in the Boltzmann 
limit the combination T 3 {dW^\2 / dV dt)^ 1 is independent of the fugacity. 
The second equality follows from the detailed balance condition, Eq. (|A.2p . 

Entropy is conserved in the expanding Universe, and in the radiation 
dominated Universe we have d(T 3 V)/dT = 0, and hence \/t s = for the 
present discussion. 1/tt describes the effect of dilution of phase space and 
can be dominant in situations where the number of the type '3' particles is 
preserved but their density decreases due to expansion of the volume. In the 
early Universe, comparing tt to the chemical relaxation time T3 provides 
the quantitative condition for freeze-out from chemical equilibrium 

T 

73 — Tr — tt- (36) 

7TI3H 

In the second equality, we have given the relevant estimate of tt for the 
early Universe [5]. 

We can now study the reaction Eq. (|26p quantitatively. In the early Uni- 
verse, the number of photons is high and is characterized by the observed 
value of 77, see below Eq. (|2U|) . and pions are also copious immediately after 
the phase transformation. These conditions, specifically chemical equilib- 
rium of the HG, mean the fugacities T f are unity and therefore the quantum 
statistics of the gases are important, leading to bose enhancement and fermi 
blocking as the case may be, recognized by Uehling and Uhlenbeck [23J. The 
importance of this effect is shown at left in Figure [JJ 

Comparing T3 to tt (l/H) at right shows that ir° remains in chemical 
equilibrium even as its thermal number density gradually decreases, consis- 
tently with falling thermal production rates. Even so at all times ttq remain 
in chemical equilibrium: this phenomenon can be attributed to the high 
population of photons, within which it remains probable to find photons 
of high enough energy to fuse into 7r°, and as the number of high energy 
photons decreases, described by photons' Planck distribution, so does the 
number of ttq which need to be maintained. 
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Fig. 7. (Color online) At left: The Boltzmann (low phase space occupancy) limit 
[dashed (blue) line] to the quantum statistical gas [solid (green) line] with fugacity 
T f — 1. At right: ir° equilibration time T3 [solid (blue) line] and Universe expansion 
time 1/H [dashed (green) line] as functions of temperature. 



One can generalize Eq. ([32]) for the three-body muon decay Eq. (j29|) and 
the 2-to-2 reactions Eqs. (f30|) , (f3T|) relevant for pion and lepton equilibra- 
tion |24] [5l H]. For example, for the three-body decay there would be three 
numerical subscripts 1,2,3 on the right, each for a distinguishable particle 
species and the product of their three fugacities appears in the equation for 
chemical relaxation time, T4, similar to Eq. (|35p . 

From this study, we learn several important characteristics of the plasma 
present in the Universe preceding and up to the time of BBN. First, tt^ 
and remain in chemical equilibrium until T ~ 6 MeV and ~4 MeV, 
respectively [5], alone due to the specific reactions considered - if and when 
other reactions are important the actual freeze-out temperature could be 
still lower. Note here that we evaluate the reaction rate for muons Eq. (|29p at 
low temperatures T <C based on the muon lifespan in vacuum 2.2010 -6 s, 
muon equilibrium is therefore preserved due to detailed balance. However, 
muon decay is a process in which helicity considerations keep the reaction 
slow, and in the dense early Universe plasma muon decay rate can increase, 
which by detailed balance argument implies that the back reaction also is 
enhances, and the chemical equilibrium is maintained longer. 

All this means that a significant number of /U is present along with 
7T° to interact with nucleons even at low (but not necessarily BBN-range 
) temperatures, as shown at left in Figure EJ A further consequence of 
the maintenance of ir, fi chemical equilibrium is a relatively large effective 
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T [MeV] 



Fig. 8. (Color online) Neutrino reaction relaxation times compared to Hubble rate. 
Solid (blue) line is for Eq. (|2"5)) while dotted (red) and dash-dotted (green) are for 
the electron- and muon-neutrino processes Eq. (|5T)) . 



degeneracy, implied also by the significant fraction of mesons in the energy 
budget, Fig. [H 

From neutrino population dynamics, we learn that relatively heavy par- 
ticles (m 3> T), in this case pions and muons, can be important influences 
on light particles' equilibration and population. At right, in Figure El one 
can see how important the pion-muon-neutrino reaction Eq. (|28p is in main- 
taining neutrino chemical equilibrium in this epoch. It appears to be more 
important than 

e + + e~ ^ v e ^ + u e>lx , (37) 

whose rate is obtained in the limit of small neutrino chemical potential 
fi u <^T [25] (shown as the dotted and dashed dotted lines). 

We have only highlighted here a few outcomes of the study to illustrate 
the necessity of studying quantitatively the chemical evolution and equili- 
bration of the early Universe plasma. More results and discussion are found 
in 0. 
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5. Conclusions 

We have shown that hadrons remain in chemical equilibrium through- 
out the evolution of the Universe. The same mechanism that allows tto to 
rapidly decay in reaction Eq. (|26p is present when two photons collide in the 
thermally equilibrated Universe, and photon 'fusion' reactions fill any miss- 
ing 7To. In order to see that there is no freeze-out of pions, a quantitative 
study is required [H [5] . 

Given that hadron chemical equilibrium is maintained, one can use meth- 
ods of hadro-chemistry to study particle abundances. We have shown how 
to compute the values of the quark and lepton chemical potentials that 
yield the observed matter-antimatter asymmetry in the Early Universe after 
hadronization and equilibration of the HG. Furthermore, we have seen that 
the non-zero chemical potentials drive charge distillation during the phase 
transformation, with the QGP and HG having negative and positive charge 
densities, respectively. Separation of baryons and antibaryons into domains 
could maintain a homogeneous zero charge density Universe, a phenomenon 
which could, e.g., play a significant role in amplifying a pre-existent, much 
smaller net baryon yield. 

We have described the chemistry of the HG in the Universe evolution pe- 
riod down to T ~ 3 MeV noting a significant yield of pions and muons which 
remain in chemical equilibrium throughout. Consequently, they remain a 
significant fraction of the particle number and energy in pre-BBN Universe. 
Their presence and reactions can impact the relic neutrino population. 

The lessons of this study are manifold, perhaps the most salient of which 
is the importance of studying the chemistry of the early Universe quark and 
hadron plasma quantitatively. Future work will be necessary to explore the 
implications of the findings here for, e.g., baryon separation, relic neutrino 
background, and thus big bang nucleosynthesis. 

Appendix A 

Fugacity Evolution Equation 



Consider the case of an unstable particle 3 with a two-body decay to par- 
ticles 1 and 2. The number of '3' particles is governed by the population 
equation 

1 dN 3 = dW 12 ^3 dWs^n , 
V dt dVdt dVdt ' [ ' ' 

which gives the change in number of particle 3 as increased by production 
(particles 1 and 2 fuse into 3) and decreased by the two-body decay 3 — > 1+2. 

In the absence of ambient populations of particles 1, 2, 3, there will be 
no production and decay will occur as in vacuum. In this case, Eq. (jA.ip 
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reduces to V~ l dN^/dt ~ — {^To) -1 N3/V , where To is the lifetime of particle 
3 in vacuum, and results in the usual exponential decay of particle 3 number. 
Here, 7 is the Lorentz boost factor relating the particle rest frame (where 
To is determined) to the observer frame; in our case, the observer is the heat 
bath of the Universe. 

The explicit forms of the rates ^y^ 3 an d d ^vdt 2 are discussed in [24"1 
[5]. Considering their structure, energy-momentum conservation and time- 
reversal symmetry of the 3^1 + 2 matrix element together imply the 
detailed balance relation 

Therefore, in chemical equilibrium, when dN^/dt = 0, Eq. (jA.ip shows that 
the fugacities must satisfy 

T 3 = TiT 2j (A.3) 

which corresponds to the Gibbs condition for the chemical potentials. 

Using Eq. (|A.2p . one can rewrite Eq. (jA.ip into the differential equation 
for T presented in Eq. (j32|) . 



Acknowledgments This work has been supported by a grant from the U.S. 
Department of Energy, DE-FG02-04ER41318; Laboratoire de Physique Theo- 
rique et Hautes Energies, LPTHE, at University Paris 6 is supported by 
CNRS as Unite Mixte de Recherche, UMR7589. 



REFERENCES 



[1] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86, 010001 (2012), see 
Section 21. 

[2] J. Letessier and J. Rafelski, Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 
18, 1 (2002). 

[3] M. J. Fromerth and J. Rafelski, "Hadronization of the quark Universe," 
|astro-ph/0211346| 

[4] I. Kuznetsova, D . Habs and J. Rafelski, Phys. Rev. D 78, 014027 (2008) 
|arXiv :0803. 15881 [hep-ph]] . 

[5] I. K uznetsova and J. Rafelski, Phys. Rev. C 82, 035203 (2010) 
|arXiv:1002.0375l [hep-th] ] . 

[6] D. J. Fixsen, E. S. Cheng, J. M. Gales, J. C. Mather, R. A. Shafer, and E. L. 
Wright, Astrophys. J. 473, 576 (1996). 

[7] A . G. Cohen, A. D e Rujula and S. L. Glashow, Astrophys. J. 495, 539 (1998) 

astro- p h/9707087. 
[8] A. D. Dolgov, Phys. Atom. Nucl. 74, 462 (2011). 



24 



Rafelski'QGU3 printed on November 20, 2012 



[9] R. Battiston [AMS Collaboration], pp741-750, in Como 2009, 11th ICATPP 
Conference on Astroparticle, particle and space physics, detectors and medical 
physics applications edited by: Claude Leroy et al, (World Scientific; Singapore 
2010). 

[10] R. Hagedorn and J. Rafelski, Phys. Lett. B 97, 136 (1980). 

[11] S . Hamieh, J. L etessier and J. Rafelski, Phys. Rev. C 62, 064901 (2000) 
|hep-ph/0006085]. 

[12] C. Bernard et al . [MILC Collaboration], Phys. Rev. D 71, 034504 (2005) 
| hep-lat/0405029 ; . 

[13] O. Philipsen, "The QCD equation of state from the lattice," larXiv:1207."5999l 
[hep-lat] . 

[14] A. Ba zavov et al. [Ho tQCD Collaboration], Phys. Rev. D 86, 034509 
(2012) |arXiv:1203.0784l [hep-lat]] ; A. Bazavov, et al. [HotQCD Collaboration] 
"Freeze-out Conditions in Heavy Ion Collisions from QCD Thermodynamics," 
larXiv: 1208.12201 [hep-lat]. 

[15] N. K. Glendenning, "Compact stars: Nuclear physics, particle physics, and 
general relativity," New York, USA: Springer (2000) 390 p 

[16] F. P. An et al. [DAYA- BAY Collaboration], Phys. Rev. Lett. 108, 171803 
(2012) |arXiv:1203. 16691 [hep-ex]]. 

[17] G. Steigman, "Neutrinos and big-bang nucleosynthesis," arXiv:1208. 0032 [hep- 
Ph]- 

[18] G. Steigman, Fortsch. Phys. 50, 562 (2002) [astro-ph/0202187| . 

[19] W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. Flannery, (Cam- 
bridge, 1992, 2ed) ISBN-0-521-043108-5. 

[20] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) 
jarXiv: 1001. 45381 [astro-ph.CO]]. 

[21] C. Greiner, P. Koch and H. Stoccker, Phys. Rev. Lett. 58, 1825 (1987). 

[22] C. Greiner and H. Stoecker, Phys. Rev. D 44, 3517 (1991). 

[23] E.A. Uehling and G. E. Uhlenbeck, Phys. Rev. 43, 552 (1933); See also 
L.P. Kadanoff and G. Baym, Quantum Statistical Mechanics, (Benjamin, New 
York, 1962). 

[24] I. Kuznetsova and J. Rafelski, Phys. Lett. B 668, 105 (2008) |arXiv:0804.3352l 
[nucl-th]]. 

[25] K. Freese, E. W. Kolb and M. S. Turner, Phys. Rev. D 27, 1689 (1983). 



